Analytical ground state for the three-band Hubbard model 
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For the calculation of charge excitations as those observed in, e.g., photo-emission spectroscopy 
or in electron-energy loss spectroscopy, a correct description of ground-state charge properties is 
essential. In strongly correlated systems like the undoped cuprates this is a highly non-trivial 
problem. In this paper we derive a non-perturbative analytical approximation for the ground state 
of the three-band Hubbard model on an infinite, half filled Cu02 plane. By comparison with 
Projector Quantum Monte Carlo calculations it is shown that the resulting expressions correctly 
describe the charge properties of the ground state. Relations to other approaches are discussed. 
The analytical ground state preserves size consistency and can be generalized for other geometries, 
while still being both easy to interpret and to evaluate. 



I. INTRODUCTION 

The-, interest devoted to the three-band Hubbard 
modeB is due to the fact that it describes the charge 
properties of a CuC>2 layer, like those found in the high- 
Tc superconducting cuprates, while still being compara- 
tively simple.0 The basic assumption which leads to the 
three-band Hubbard model is that the only relevant or- 



bitals are the Cu 3cL 



and the O 2p x and 2p y 



bitals. This granted, the CuC>2 layer may be described 
by a lattice with one Cu and two O sites per unit cell with 
hybridization between nearest neighbor Cu-0 pairs and 
0-0 pairs. In the hole picture the three-band Hubbard 
Hamiltonian reads 



H — Ho + Hi , 

H = Aj2n p ja + U d J2nt r ni l . 

ja i 
Hi = t pd ^ ftpd(p]a d ™ + h - C -) 

+ t pp 4 ! ppP\ a Pj'(J ' 



(la) 
(lb) 



(lc) 



where d\ a (p] a ) create a hole with spin a in the i-th 
Cu 3d orbital (j-th O 2p orbital), while nf a (n^ a ) are 
the corresponding number operators. Hq is the atomic 
part of the Hamiltonian with the charge-transfer energy 
A and the on-site Coulomb repulsion Ud between Cu 3d 
holes. Hi represents the hybridization of Cu 3d and O 
2p or bitals (hopping strength t p d) and of O 2p orbitals 
(hopping strength t pp ). The factors <#Jj and (j^J p give the 
correct sign for the hopping processes^, and (ij) denotes 
the summation over nearest neighbor pairs. 

In Eq. ( JTa| ) only the most important Coulomb repulsion 
Ud is included, while O on-site and inter-site Coulomb 
repulsions have been neglected. For explicit calculations 
the following typical set of values for the parameters in- 
volved in Eq. (Qq) will be usedel 



Strong correlations due to Ud are the reason why ground 
state properties of Hamiltonian H can be calculated 
only approximately and/or on finite clusters. Besides 
analytical approaches like, e.g., different approxima- 
tions for dynamical Green's functionscl mostly numerical 
simulationsO have been applied to the three-band Hub- 
bard model. The aim of the present work is to derive an 
analytical approximation for the ground state of an infi- 
nite system at half-filling (i.e. one hole per Cu site) which 
correctly describes charge properties and is still compara- 
tively easy to evaluate. The resulting approximation does 
not only allow for a calculation of ground-state proper- 
ties of the three-band model (|la|). It also provides a 
framework for the investigation of excitations. Further- 
more, the approach is sufficiently general to be applied 
not only to a Cu02 plane but also to different geometries 
like, e.g., that of a CUO3 corner-sharing chain.El 

Starting point of the approximation is a Neel-ordered 
ground state of the atomic Hamiltonian Hq, which is de- 
noted by IV'o)- Due to fluctuations induced by Hi this 
atomic ground state \ipo) differs from the full ground state 
|\I>) of Hamiltonian H. In Sec. Ill it will be shown that a 
perturbative treatment of these fluctuations breaks down 
for parameter values which are in the physically relevant 
range. Therefore, ground-state fluctuations have to be 
treated in a non-perturbative way. In the following we 
will present a systematic and non-perturbative scheme to 
introduce these fluctuations on the background of |^>o}- 

The paper is organized as follows. In Sec. II the general 
formalism is presented. As an illustration in Sec. Ill this 
formalism is applied to the (exactly solvable) problem 
of a single CUO4 plaquette. The approximative ground 
state of an infinite, half-filled Cu02 plane is developed in 
Sec. IV, and ground-state expectation values are evalu- 
ated in Sec. V. In Sec. VI the results of the analytical ap- 
proach are compared to Projector Quantum Monte Carlo 
simulations. The conclusions are presented in Sec. VII. 
Finally, a more detailed justification of the approach to- 
gether with a discussion of its relationship to the cumu- 
lant formalism!] is given in the Appendix. 



A = 3.5 eV, U d = 8.8 eV, 
tpd = 1.3 eV, tpp = 0.65 eV. 



(2) 
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II. GENERAL FORMALISM 

The basic idea is to start with a state \ipo) which is 
a first approximation to the full ground state IvP). In 
the present case we will choose | ipo } to be a Nccl-ordcrcd 
ground state of the atomic Hamiltonian Hq, Eq. (|Tb|). 
In |^o) every Cu site is singly occupied and all O sites 
are empty. Next, fluctuations on the background of j^o) 
are introduced. These ground-state fluctuations are de- 
scribed by fluctuation operators F a (introduced below in 
more detail) which approximately transform state \ipo) 
into the full ground state \^f). Under quite general con- 
ditions it can be shownB that the transformation leading 
from ) to \^) has to be of exponential form 



|*) = exp 




(3) 



The parameters A Q are fluctuation strengths of the fluc- 
tuation operators F a . They are determined using the set 
of equations 



0=(*\[H,F*] I*) ,a = l,2,. 



(4) 



Eq. (Q) follows from the condition that is an eigen- 
state of the full Hamiltonian H, Eq. (|Ta|) . From Eq. (||) 
all ground-state properties can be evaluated using 



(A) 



(5) 



The ground-state energy, for instance, is calculated from 
Eq.(g) with A = H. Equations (||) and (||), together with 
an appropriate choice of fluctuation operators F a , con- 
stitute the formal framework to be used in the remainder 
of this work. The above formalism will be applied to an 
exactly solvable problem in Sec. III. This serves both as 
an illustration of the approach, and as an indication to 
which extent a perturbative treatment of ground-state 
fluctuations is possible. 

We finally remark that Eqs. (|3|) and (Q) are closely re- 
lated to the cumulant formalism (see Appendix) . There- 
fore the approach presented in this work preserves 
size consistency. Thus, for instance, the approximated 
ground-state energy remains an extensive quantity. 



III. APPLICATION TO A SINGLE PLAQUETTE 

As an example, the method presented in the last sec- 
tion is now used to find an approximate expression for 
the ground state of Hamiltonian ( |Ta| ) for the case of a 
single CuC>4 plaquette occupied by a single hole. In this 
case the unperturbed ground state \tpo) of H is a state 
in which the Cu site is singly occupied while the four O 
sites are empty. For reasons of symmetry we may use 
a single fluctuation operator F± which describes fluctua- 
tions of the hole from the Cu site to the four surrounding 
O sites j, i.e. 



(6) 



j a 



where (j> pd are the phase factors introduced in Eq. (Jlc 
According to Eq. (||) the full ground state of a hole on a 
single plaquette is expressed by 



|*) = exp(AiFi) \ip ) 
The norm of this state is 

(* | *) = 1 + 4Xj . 



(7) 



(8) 



The fluctuation strength Ai in Eq. (Q) is determined from 
condition (||) 



o = (*l \h,f}] I*) . 



(9) 



Non-vanishing contributions in Eq. (Q) arise only from 
terms up to order A^. The following quadratic equation 
for Ai is obtained 



= U pd - 4 (A - 2t pp ) Ai - m pd \\ 



(10) 



When the positive solution for Ai is used in Eq. (^) one 
obtains the exact ground state. The ground-state energy 
Eq is calculated using Eq. (j|) 



Eg — — 4£ P dAi 
1 " 
~ 2 



2tr, 



(A - 2t p 



(4t pd ) 



(11) 



For parameter set @ a value of Ai = 0.33 results. Notice 
that Eq. ( p"T| ) contains only a reduced effective charge- 
transfer energy A — 2t pp . The Cu-occupation number 
(n 0u ) = (*l n d |*) (* | is given by 



("-Cu) 



(V>o| exp ( AiF/ ) n d exp (Ai-Fi) |^o> 



(*l*> 

(Vto| (l + AiF+) n d (1 + A1F1) |Vo) 



(*l*> 



1 



l + 4Af 



(12) 



From this result one may conclude that a perturbative 
treatment of F\ fluctuations (i.e. an expansion in Ai, see 
Appendix) , also in infinite systems, is in general not pos- 
sible. Typically, Ai is of the order 1/2. An expansion of 
Eq. ( p^| ) in Ai, however, diverges for Ai > 0.5. Condition 
Ai = 0.5 is equivalent to a vanishing effective charge- 
transfer energy A — 2t pp = 0, and to a Cu-occupation 
number of 1/2. At this point state |i/ , o) ceases to be a 
good approximation of the exact ground state |\&). This 
divergence has been observed previouslya, although its 
origin was unclear at that time. 
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IV. APPLICATION TO AN INFINITE CUO a 
PLANE 

We now apply the formalism presented in Sec. II to the 
geometry of an infinite CuC>2 plane. State j^o) is again 
the Neel-ordered ground state of the atomic Hamiltonian 
.Ho, Eq. (lb). Let us introduce appropriate fluctuation 
operators F a . First, operator F\ from Eq. (^) is general- 
ized to all N Cu sites i 



F, 



E 



where the sum is over the four O sites j which surround 
Cu site i. The remaining operators are constructed in 
accordance with the following principles: (i) All opera- 
tors describe derealizations of a hole initially located at 
Cu site i. (ii) The final site in the process is reached 
via the shortest path accessible by Cu-0 hopping pro- 
cesses, (iii) A summation over all equivalent final sites 
is taken, (iv) The signs of the hopping processes are 
chosen to be the negative of the phases 4> l ^ d in Hamilto- 
nian ([lc]) (which guarantees non- negative values for the 
fluctuation strengths X a ). 
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FIG. 1. Final sites of the fluctuation operators Fi tC t- Cu 
and O sites are symbolized by squares and circles, respec- 
tively. Arrows show spin orientation and position of the 
holes in the atomic ground state. The fluctuations Fi lCl start 
from Cu site i and lead to the final sites sites labelled by 
a = 1, ...4. The area shown is a quarter of the full area 
accessible to the fluctuations. 

Figure 1 shows final sites reached by fluctuation oper- 
ators Fi t0l , a — 1, . . .4. For reasons of symmetry only a 
quarter of the allowed fluctuation range is shown. Fluctu- 
ation Fi r 2, for instance, describes the hopping of the hole 
from Cu site i, via O site j, to the four nearest-neighbor 
Cu sites k 



F 



jka 



1 



4;A> 



Final states with singly or doubly occupied Cu sites differ 
by the Coulomb energy Ud- Since Ud is large we have to 
distinguish between these two cases. Therefore we split 
Fi t 2 into two operators which describe a process leading 
to a singly or doubly occupied Cu site, respectively 



i^ = -£(l-ny (l-n%)dld i<T , 

jka 
jka 

Note that it is not necessary to introduce a fluctuation 
operator which leads to the nearest neighbor Cu sites in 
diagonal direction (e.g. the Cu site without a label in 
Fig. 1). Due to the Pauli principle fluctuations to these 
sites are largely excluded because of antiferromagnetic or- 
der. The neglect of fluctuations leading beyond the range 
shown in Fig. 1 will be justified a posteriori. It will be 
shown that the fluctuation strengths X a decrease rapidly 
with increasing length of the fluctuation processes. 

According to Eq. (||) the ground state has the following 
form 



I*} = exp 




(13) 



where a denotes the 5 fluctuation operators described 
above. Because of translational symmetry the parame- 
ters X a do not depend on the Cu-site index i. To sim- 
plify Eq. (|l^) we approximately factorize the exponential 
function with respect to the fluctuations Fi \ 



|*) = exp ^2 X aF ia exp 

\ia>l / \ i> J 

This approximation amounts to the assumption that far- 
reaching fluctuations i^ ia >i occur on the background of 
-F^i-fluctuations which in turn are influenced only indi- 
rectly (i.e. via Aj) by the former. The second exponential 
function in Eq. (fhfl) exactly factorizes with respect to i' 



|tf ) = exp (J2 A « F » II + Al ^M) l^>) ■ ( 15 ) 

\IQ>1 / V 

In Eq. (|l5|) every hole may fluctuate over a total range 
of five plaquettes each. Notice that all holes fluctuate 
simultaneously. This leads to a multitude of many-body 
effects, i.e. the fluctuation of a hole depends on the con- 
figuration of other holes. Basically there are three types 
of many body effects which are exemplified in Fig. 2. 
First, due to the Pauli principle the fluctuation of a hole 
may be blocked by the presence of other holes with the 
same spin, as in the fluctuation process labelled (a) in 
Fig. 2. Second, there are processes in which holes with 
the same spin change place, cf. process (b). In the follow- 
ing we will call these processes site-changing processes. 
Third, there are strong correlations due to the Hubbard 
Ud on doubly occupied Cu sites, as in process (c). One 
common feature of all these many-body effects is that 
they suppress fluctuations. 
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FIG. 2. Examples for many body effects. There are three 
types of effects: (a) processes which are excluded by the 
Pauli-principle, (b) site-changing processes, and (c) correla- 
tions due to the Hubbard Ud on doubly occupied Cu sites. 
The approach presented in this work accounts for all effects 
shown here. 

This multitude of many-body effects makes an exact 
evaluation of expectation values using Eq. ( |l5| ) impossi- 
ble. Further approximations are therefore necessary. Let 
us consider processes in which two or more holes simul- 
taneously leave their original plaquette. The fluctuation 
strengths X a for such far-reaching fluctuations turn out 
to be small compared to Ai . Therefore it should be pos- 
sible to neglect the many-body effects arising in these 
processes (except for site-changing processes in diagonal 
direction, see below). In the case of F^d fluctuations, 
for example, we neglect the possibility that the O site 
j between the Cu starting and final sites i and k may 
already be occupied by a hole with the same spin. This 
amounts to the following simplification 



Fa 



2d 



jka 



In this way all of the aforementioned processes are in- 
cluded, some of them however only in a simplified way 
(i.e. by neglecting many-body effects). 

In addition to all many-body effects which are due 
to processes where only one hole leaves its original pla- 
quette, we furthermore take account of all site-changing 
processes in diagonal direction, see Fig. 2(b). The sup- 
pression of charge fluctuations due to the diagonal sites 
turns out to be of great importance. On the other hand, 
site-changing processes involving next-nearest Cu neigh- 
bors in horizontal or vertical direction can be neglected 
for the following reason: During these processes the paths 
of the holes have to cross at the intermediate O sites (i.e. 
sites 1 and 4 in Fig. 1). However, since the holes have 
the same spin, they have to avoid each other due to the 
Pauli principle. Site-changing processes in horizontal or 
vertical direction are therefore unlikely. 



V. EVALUATION OF EXPECTATION VALUES 

We now evaluate expectation values with state Jl5|). 
Using the above approximations the norm of this state is 

<*|*> = 

generalization of Eq. (H) 



where N is the number of Cu sites and, in 



^2 z a p a X 2 a 



(16) 



z a is the number of equivalent final sites of the given 
process (e.g. zi s = 4), and p a is the probability that the 
configuration of the other holes makes the process possi- 
ble. This probability is defined by 



Pa 



(tt|^,q|tt) 
(* | <F) 



(17) 



where Pi ia is a projection operator on all configurations 
which allow for process a. For example, Pi.2 S projects 
on states in which the target Cu site of fluctuation Fi t 2 S 
is empty, whereas P\ H i& is the projector on states with a 
singly occupied final site. Due to translational symme- 
try the probabilities p a do not depend on the site index 
i. Obviously, Eqs. ( Jilf ) and (|l7f ) have to be solved self- 
consistently since \^) in Eq.(jl7|) depends on the param- 
eters A Q . In the case of p2 S , for instance, we obtain by 
explicit calculation 



P2s 



{MP ias m 

V N 

, ZaPaX a 



V 



= 1- l/v . 
In an analogous way one finds 

Pi = 1 , P2s = 1 ~ 1/f , P2d = I/V , 
P3 = 1 - 2\l/is , p 4 = 1 - A?/l/ • 



(18) 



The interpretation of Eq. ( |i8| ) is straightforward, pi = 1 
holds since we assume that far-reaching fluctuations oc- 
cur on the background of i^.i-fluctuations. l/v is the 
probability to find a given hole at its original Cu site. p2d 
is therefore the probability that a target Cu site is singly 
occupied. This is a necessary prerequisite for the fluc- 
tuation process Fid which leads to a double occupancy, 
cf. Fig. 2(c). p2s, on the other hand, is the probability 
that a target Cu site is empty, as required for fluctuation 
process i 7 ^. The probability to find a given hole at a 
specific O site on its original plaquette is \\/v. Thus p± 
is the probability that the target O site of fluctuation F4 
is not blocked by the hole of same spin which resides on 
the neighboring Cu site, cf. Fig. 2(a). Analogously, p$ is 
the probability that the target O site of fluctuation F 3 is 
not blocked. The additional factor 2 in p% (as compared 
to P4) is due to site-changing processes, cf. Fig. 2(b). 
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The fluctuation strengths X a are calculated using 
Eq. (0) for an arbitrary site i = 



= (*l 



I*) 



(19) 



One obtains the following nonlinear system of equations 

= {Eg — A + 2t pp ) Ai + t p d + t p dMsP2s 

+t p d^2dP2d + 2tp P \ 3 p 3 - 2t P p\\\- i /v , (20) 
= E G X 2s + -it pd X 2 dXi / v + t pd \iP2s 

+2tpdMP3P2s + tpd,XiPiP2s , (21) 

= {2E G - U dd ) \2d + ^ pd \ 2s X 1 

+t p d^i + 2tp d \zPz + t p d\±Pi , (22) 
= {Eq - A + t P p) X 3 p 3 + tppXips 

+tpdX2sP2sP3 + t pd X 2 dP2dP 3 

+t pp X A p 3 - {t pd + 2t P p\i) \\\ 3 /v , (23) 

= [Eg - A) A 4 P4 + tpd^2sP2sP4 + t pd \2dP2dP4 

+2t pp X 3 p 3 - {tpd + 2tp P X 1 ) AiA 4 /^ , (24) 

where Eg — —4t p d\i is the ground-state energy per Cu 
site (see Eq. (pq)). This system of equations, together 
with Eqs. (16|)and Jlq ) can be solved self-consistently 
for all A Q , p a and for v. The solution with the lowest 
value of Eg is then used in Eq. (|l5|). In the case of 
A Q = for all a > 1 Eq. @ reduces to Eq. @ for the 
single plaquette. 
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FIG. 3. Derealization probability as a function of fluctu- 
ation length. The graph shows the derealization probability 
of a hole which originates from Cu site i. The probability is 
summed over the sites displayed beneath the bars, and over 
all equivalent final sites. The hole remains predominantly 
on its original plaquette. Derealization beyond the near- 
est-neighbor plaquette is negligible small. The probability 
has been calculated using Eq. ( |is| ) for parameter set (Q). 

Figure 3 shows the delocalization probability p a )? a /v 
of a given hole summed over equivalent final sites as a 



function of fluctuation length for parameter set (g). In 
order to demonstrate the convergence of the results addi- 
tional fluctuations have been introduced which lead be- 
yond the fluctuation range shown in Fig. 1. The contribu- 
tion of these additional fluctuations to the ground-state 
energy amounts to less than 0.1 percent. No significant 
delocalization beyond the nearest neighbor plaquette oc- 
curs. A similar observation is made when other model- 
parameter values are chosen within the range which is 
relevant for cuprate compounds. These results retrospec- 
tively justify the neglect of far-reaching fluctuations and 
many-body effects. Notice, however, that the neglect 
of many-body effects allows for unphysical fluctuations 
which may decrease the calculated ground-state energy 
below the exact value. Thus, in contrast to an exact 
evaluation of Eq. (^|), our approximate solution does not 
guarantee an upper limit to the exact ground-state en- 
ergy. 

From ground state ( |15|) all expectation values are eas- 
ily evaluated using Eq. (||) . The ground-state energy per 
Cu site, occupation numbers, and double occupancies of 
Cu and O sites are 



E G 

("-Cu) 

{dcu) 
{no) 
(do) 



-AtpdXi , 



- (l + 4A2 s p 2a + 4A£ I p M ) , 
-(4ALpm) , 
£ (A? + 2A|p 3 + X 2 4 pa) , 
j (no) ■ 



(25) 
(26) 

(27) 

(28) 

(29) 



The number of holes is conserved, i.e. (ncu)+2 (no) = 1. 
By comparison with Quantum Monte Carlo calculations 
it will be shown in the next section that these results 
correctly reproduce the charge properties of the ground 
state. However, magnetic properties like the reduction 
of sublattice magnetization due to fluctuations are only 
partly described. We have neglected many-body effects 
in processes in which two holes simultaneously leave their 
original plaquette. Therefore, no spin-flip effects are in- 
cluded. Thus, in the (Heisenberg-) limit of infinitely large 
A and Ud ground state ( |l5| ) reduces to the Neel state. 



VI. DISCUSSION OF THE RESULTS AND 
COMPARISON TO QUANTUM-MONTE CARLO 
SIMULATIONS 

We have carried out numerical simulations of the three- 
band Hubbard model ([la]), using |the Projector Quantum 
Monte Carlo (PQMC) algorithn£3, in order to compare 
them with the analytical result, Eq. (15). In the PQMC 
approach the ground state of a finite cluster is projected 
out from a suitable trial state by applying the exponential 



operator e ^ onto the trial state in the limit j3 — > oo. 
However, in numerical calculations only finite values of 
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the parameter /3 are accessible. Therefore one has to 
check convergence of the results with respect to (3. Fur- 
thermore one has to account for possible finite size effects. 
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FIG. 4. Convergence of the PQMC calculations with re- 
spect to /3 and the system size. The ground-state energy Eg 
per Cu site for clusters of 4 x 4 and 6x6 plaquettes is shown 
as a function of inverse linear system size L. The solid and 
broken error bars are the results for (5 = 4 and /3 — 8, respec- 
tively. The parameters are those of set (^). 

In the present study two different cluster sizes have 
been investigated. Firstly, we have used a system con- 
sisting of 4 x 4 plaquettes (i.e. 48 sites). This is the 
smallest cluster which allows for periodic boundary con- 
ditions while still being fully two-dimensional. Secondly, 
we have studied the next largest cluster, a system of 6 x 6 
plaquettes (i.e. 108 sites). We have used a mean-field ver- 
sion of the analytical ground state (|l5|) as trial state. No 
sign problem occurred. 
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FIG. 5. Comparison of analytically calculated ground-state 
energies and Cu-occupation numbers (solid lines) with the re- 
sults of PQMC simulations for a cluster of 4 x 4 plaquettes 
(error bars connected by broken lines). As a function of A 
the plots show the ground-state energy Eq per Cu site (upper 
graph), and the Cu-occupation number (ncu) (lower graph). 
The other parameters are those of set (g). 



It turns out that the results obtained for the 4x4 sys- 
tem with j3 = 4 are already reasonably well converged 
with respect to both f3 and system size. As shown in 
Fig. 4 for the case of the ground-state energy per Cu site 
the error bars for different values of [3 overlap and the 
values of the 4x4 system differ only slightly from those 
of the 6x6 system. For this reason we restrict our sim- 
ulations to a system of 4 x 4 plaquettes with /3 = 4 and 
compare the results with the analytical approach. 

In Figs. 5 and 6 the ground-state energies Eq per 
Cu site and several occupation numbers calculated us- 
ing Eqs. (p5|)-(p9|) are compared to the results of the 
PQMC simulations. The number of holes is conserved in 
both approaches. Thus the O-occupation number (no) is 
a function of the Cu-occupation number (nc u ), and the 
former is therefore not shown. While the values of the 
model parameters are those of set (|j), the charge-transfer 
energy A is varied covering the range of very large charge 
fluctuations (A = 1.5 eV) to fairly small charge fluctua- 
tions (A — 4.5 eV). In general there is a good agreement 
between the analytical and numerical results, especially 
for larger values of A. 
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FIG. 6. Comparison of analytically calculated double oc- 
cupancies (solid lines) with the results of PQMC simulations 
for a cluster of 4 x 4 plaquettes (error bars connected by bro- 
ken lines). As a function of A the plots show the Cu-double 
occupancy (dcu) (upper graph), and the O-double occupancy 
(do) (lower graph). The other parameters are those of set 
(I)- 

With increasing A both Eq and (ncu) increase while 
the O-double occupancy (do) decreases. This behavior is 
due to the suppression of fluctuations for larger values of 
the charge-transfer energy. For smaller A the analytical 
value for the ground-state energy lies below the PQMC 
result. This can be explained by the neglect of many- 
body effects in Eq. ([l5]) which allows for more unphysical 
fluctuations when A becomes smaller. These fluctuations 
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also lead to values for (ncu) which are slightly larger than 
the PQMC result. However, even for A = 1.5 eV the rel- 
ative deviation for both Eg and («c u ) amounts to less 
than 3%. 

For larger values of A the decrease of (do) with in- 
creasing A is about two times as large as the decrease of 
(dcu)- The reason for this weaker dependence of (dcu) 
on A is that an increase in the charge-transfer energy af- 
fects (dcu) only indirectly by reducing the effective Cu- 
Cu hopping, while - in contrast to (do) - the on-site en- 
ergy of the final site is not changed. Both analytical and 
numerical results show a maximum in the Cu-double oc- 
cupancy (dcu) for intermediate values of A. This may be 
interpreted as the point where A is already sufficiently 
large to force holes from O sites onto already occupied Cu 
sites but still not large enough to suppress the effective 
Cu-Cu hopping. 



linear combination of different factorizations of expecta- 
tion values. For example, for two operators Ai and Ai 



(<p\AiA 2 \ X ) c 



( V >\A l A 2 \ X ) (^|^i|x)(vl^ 2 |x) 



(<p I xY 



Here and in the following we always assume that the 
states involved in a cumulant have non-vanishing over- 
lap, i.e. (tp | x) 0- One of the attractive features of 
cumulants is that they preserve size consistency.!] 

The cumulant formalism for the calculation of ground- 
state properties may be formulated as follows. By appli- 
cation of an operator Q within the cumulant ordering an 
approximate ground state |^ ) can be mapped onto the 
full ground state |\&) 



|*) c = f! |^ } c = cxp ( ^jA a F Q ) 



(Al) 



VII. CONCLUSION 

Summing up, we have derived an analytical approx- 
imation, Eq. (|l5|), for the ground state of the three- 
band Hubbard model (|la|) on an infinite, half filled Cu02 
plane. The approach uses fluctuation operators F a and 
fluctuation strengths A Q which have a clear physical inter- 
pretation. The parameters contained in Eq. ( |l5| ) are de- 
termined self-consistently by solving a nonlinear system 
of equations. While the approach is non-perturbative 
and conserves size consistency, expectation values with 
the approximate ground state are still easy to evaluate. 
By comparison with Projector Quantum Monte Carlo 
simulations we have demonstrated that Eq. (|l^) gives 
a reliable description of charge properties covering the 
range from small to very large charge fluctuations. Equa- 
tion (|l5| ) can be generalized for other geometries. Fur- 
thermore, due to the use of fluctuation operators our ap- 
proach provides a natural framework for the calculation 
of charge excitations, for example by using projection 
technique. This will be demonstrated in a forthcoming 
publication^] 
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APPENDIX: RELATION TO THE CUMULANT 
FORMALISM 

The approach presented above can be formulated in 
the framework of the pimulant formalism!] The cumu- 
lant expectation valudiil of a product of operators is a 



The parameters, X a are determined using the following 
set of equations! 



0=(<p\FiHn\iP o ) c , 



(A2) 



for alLa and for an arbitrary state \tp). These equations 
followE£l from the condition that f2 \ipo) c is an eigenstate 



of H. The exponential function in Eq. (Al) should be un- 
derstood in terms of a series expansion in which the op- 
erators F a are subjected to the cumulant ordering. From 
Eq. (Al) ground state properties can be calculated using 



(A) = (<p\ An l^o) 



(A3) 



Next we show that Eq. (j^) and Eq. (Q) can be derived, 
from the above equations when the following identity££l 
is used 



( V \e F] Ae F \ X ) c = {e F <p\A\e F x) C 



(A4) 



Equation (A4) holds for all operators F and A. It allows 



to remove the exponential functions from the cumulant 
orde ring and apply the m di rectly onto the states. Using 
Eq. (|A4j), Eqs. (ph and (Al) can be directly transformed 
into each other. In Eq. (A2), on the other hand, we 
choose \<p) = O \ipo) and use Eq. (A4) to obtain 



0= (Q^o\FlH\n^ ) c ■ 

If |n^>o) is an eigenstate of H this equation is equivalent 
toEq. |). 

We conclude by pointing out some of the advantages of 
Eqs. ([|) and (^) as compared to other possible approaches 
within the framework of the cumulant formalism. First, 
the fact that the exponential operator Q has been trans- 
ferred onto the state \ipo) amounts to a summation of all 
orders in the A Q . Thus, our approach is non-perturbative 
and avoids (possibly divergent) series expansions of ex- 
pressions like Eq. (|l2|). The divergence of these series, 
for Ai > 0.5 is equivalent to the violation of conditionE£l 
\(ipo \^}\ 2 > 1/2- Consequently, when Ai > 0.5, within 
the cumulant ordering no operator fl exists for state |^>o)- 
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These difficulties do not occur when fl has been removed 
from the cumulant ordering. Second, since in Eq.(^|) the 
full ground state |^} appears as both bra and ket vector 
more flu ctuations are taken into account. If, for example, 
Eq. (A2) with \(p) — \ipo) is used instead of Eq. (|j), one 
always obtains a vanishing value for A2 S - 
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